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Abstract. We study the following quasilinear elliptic system for all i = 1, • ■ ■ , m 

-div(#'{|Vui|2)Vui) =in R" 

where u = : R" — t R"* and the nonlinearity Hi(u) g C^(R'") — >■ R is a general nonlinearity. Several 

celebrated operators such as the prescribed mean curvature, the Laplacian and the p-Laplacian operators 
fit in the above form, for appropriate $. We establish a Hamiltonian identity of the following form for all 
a;„ e R 


/ 

v/R 


R ^-1 




<i.(|Vui|2) -#'(|vui|2) 


- H{u) dx' = C, 


where x = (x\xn) € and H is the antiderivative oi H = This can be seen as a counterpart of 

celebrated pointwise inequalities provided by Caffarelli, Garofalo and Segala in [13 and by Modica in [^ . 
For the case of system of equations, that is when m > 2, we show that as long as 

(2s^'(s)' 


a > a* := inf < - 
“ s>Q { 


$(s) 


the function •= .J-a Sb — is monotone nondecreasing in r. This in particular 

implies that for the prescribed mean curvature, the Laplacian, the p-Laplacian and operators the function 
is monotone when q > a* = 2, q > a* = 2 and o; > q* = p, respectively. We call this a weak 
monotonicity formula since for m = 1 it is shown in [l^ that /a(r) is monotone when a > 1, under certain 
conditions on 

We prove De Giorgi type results for 7f-monotone and stable solutions in two and three dimensions when 
the system is symmetric. The remarkable point is that gradients of all components of solutions are parallel 
and the angle between vectors Vu-i and Vuj is precisely arccos ^ ^ ^ ^. In addition, we provide an optimal 

Liouville theorem regarding radial stable solutions of the above system with a general nonlinearity when the 
system is symmetric. We announce several natural open problems in this context as well. 
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1. Introduction 


In [17|, Caffarelli, Garofalo and Segala studied the following class of quasilinear equations arising in 
geometry 

(1.1) div($'(|Vu|^)VM) = f{u) in K", 


where / is a and $ G C'^(K+) satisfies certain conditions that follow. Note that for the case of 

d>(s) = s the above equation is the standard semilinear elliptic equation. The prescribed mean curvature 
equation and the p—Laplacian equation, i.e. 


( 1 . 2 ) 

div 

(1.3) 

div 


Vu 


Vl + |Vu|^ 


= fiu) in 


div (e+lVu^) —Vu =/(u) in K" for e > 0 , 


respectively, fit in the form of (HID where $ is given respectively by 

(1.4) $(s) = 2(yrT7-l), 

(1.5) $(s) = ^ (^(e + s)5 - . 

Throughout this paper we shall assume that $(s),$'(s) and 4>'(s) + 2$"(s)s are positive when s > 0. 
In addition, without loss of generality let 4>(0) = 0. Borrowing notations from 0 , we shall refer to the 
following conditions often in this paper. Note that $ in (11.41) and (11.51) satisfies these conditions, respectively. 
Condition (A). There exist positive constants Ci,C 2 and e > 0 such that $ G C'^(R+) and for every 
r;, C G K" 

(1.6) Ci(e+M)-i<<i>'(hn<C2(e+M)-\ 

m 

(1-7) ^1(6+ |77|)"^|CT < XI < <^2(6+ b|)”^|CT> 

i,J = l 

where C,' = (Ci Cn+i) ^ is orthogonal to the vector (— 77 , 1 ) G 


Condition (B). There exist p > 1, e > 0 and positive constants Ci,C 2 such that $ G C'^(K+) that 
for every ? 7 , C € M", 

(1.8) Ci(e + < ^^(e + |r7l)"-^ 

m 

(1-9) Ci(e+ IryD^^-^ICp < X < ^^(e + 

where aij in ISTM and dni) are given by 




IJ . 


( 1 . 10 ) 

One of the main results provided in 0 is the following pointwise inequality. Note that this is a counterpart 
of the pointwise estimate given by Modica in ^ for the case of $(s) = s. 


Theorem A. Suppose that f G C^(R) with F > 0 and suppose that one of the following conditions hold 

(i) Condition (A) holds and u G Wdoc is a solution to U.l\} 

(ii) Condition (B) holds and u G C'^(R"') (~lL°°(R") is a solution to U.l\) and in additon |Vu| G L°°(R"). 
Then for every x 


( 1 . 11 ) 


2$'(|Vun|Vup - 4>(|Vun < 2F{u). 


In particular, the following pointwise estimates hold for specific $. 
• Suppose that $(s) = s, then 


( 1 . 12 ) 


|Vup < 2F{u) in R”, 


where u is a bounded solution of the semilinear equation Au = /(it) in R", provided by Modica in 

ii]. 
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(1.13) 


Let $(s) = 2{y/Tl^ — 1). Then 


V^l + |VnP- 1 

x/l + |VuP 


< F{u) in 


for bounded solutions of the prescribed mean curvature equation that is div 


I 

I v'l+iv^p 


= /(w) in 


• Suppose that <i)(s) = ^S 2 . Then 


(1.14) 


|Vu|P < 


p- 1 


F{u) 


where u is a bounded solution of the p-Laplace equation div (iVul^’ ^Vu) = f{u) in 
We study classical solutions of the following quasilinear system of equations 

(1.15) -diy{<i>'{\Vu,\^)Vui) = H,{u) inR", 

where u = (m*)™i : M" ^ R"* and Hi{u) G (^^(R™) ^ R for all i = 1, 
variational structure and the associated energy functional is given by 


m. The above system has 


(1.16) 


r 1 


where H is defined such that diF[{u) = Hi{u). Throughout this paper we use the notation u = 

F[{u) = {F[i{u))^i and djFli{u) = ■ We assume that diF[j{u)djHi{u) > 0 for 1 < i < j < to. The 


next definition is the notion of the symmetric systems, introduced by the author in 3^. Symmetric systems 
play a fundamental role throughout this paper when we deal with the energy functional given in (I1.16|) and 
when we study system (11.151) with a general nonlinearity H{u). Note that for the scalar equation case, that 
is when to = 1, (11.151) is clearly symmetric. 


Definition 1.1. We call system hl.l5\] symmetric if the matrix of partial derivatives of all components of 
H given by 


(1.17) 

is symmetric. 


:= (9.77, (n))™.^!. 


Hamiltonian identities are quite well-known in both mathematics and physics as important tools to study 
qualitative behaviour of entire solutions of differential equations and systems. They often directly or indi¬ 
rectly lead to certain properties which could be of great importance in the fields as well, such as monotonicity 
formulae. Consider the following symmetric system of ordinary differential equations that is a particular 
case of (11.151) . 

(1.18) —u'f = diH{u) in R. 

It is straightforward to see that the following Hamiltonian identity holds for solutions of (11.181) 


(1.19) 


2 ^ 

i=l 


u'f + H{u) = C in 


where C is a constant. Equivalently, one can rewrite (11.181) in the form of a first order Hamiltonian system 

-u[ = dy^H{u,v) in R, 

—n' = —duiH{u,v) in R, 

where H{u,v) = ^ + H{u). Note that H{u,v) = C on trajectories of solutions. These equations 

generalize Newton’s third law that is F = ma to system of equations where the momentum is not simply 
mass times velocity. The Hamiltonian H{u,v) normally represents the total energy of the system. We refer 
interested readers to [3, for some original information regarding physical meaning of the system and to 
MEl and references therein for variational theory of Hamiltonian systems. 
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Gui in [39j considered the gradient system —Aui = diH{u), that is a higher-dimensional counterpart of 
(11.181) . and established the following Hamiltonian identity, 


( 1 . 20 ) 




IIL 

-T 
2 ^ 

i=l 




dx' = C, 


for X = {x',Xn) G X R. In this paper, we provide an extension of this inequality for solutions of 

quasilinear symmetric system (jl.l5|) . One might expect, at the first glance, that just replacing derivative 
terms with could simply give the Hamiltonian identity for solutions of (11.151) . 

Instead, the identity follows the structure of the pointwise estimate provided by Caffarelli, Garofalo and 
Segala in 0 and it is of the form 


( 1 . 21 ) 


E 


-$(|Vu,|2) -$'(|Vu,|2) \dx^u,\ 


- H{u{x)) dx' = C. 


Note that when $(s) = s, the identity (11.211) recovers (ll.20p . This then explains why the difference of partial 
derivatives that is \'Vx'Ui\'^ — \dx„Ui\'^ appears in (ll.20p . 

If we set <I> to be the ones given in (d and m we can have the Hamiltonian identity for the prescribed 
mean curvature equation and the p-Laplacian equation, respectively. Let us mention this remarkable point 
again that the Hamiltonian identity ()1.21l) has a very similar structure as pointwise estimates ()l.lll) and 
(11.121) . provided by Caffarelli, Garofalo and Segala in 171 and by Modica in [i^. Therefore, (11.211) can be 
seen as a counterpart of (11.111) for the case of system of equations, i.e. m > 2. 

The Hamiltonian identity (11.211) motivates us to look for a monotonicity formula for solutions of (11.151) . 
So, set 


( 1 . 22 ) 


4(r) := ^ [ 
^ JB, 


^$(|vu.n-2H(w) 


For the case of scalar equations, that is when m = 1, it is proved by Caffarelli, Garofalo and Segala in [171 
that when $ satisfies one of conditions (A) or (B) then the function Ia{r) is monotone nondecreasing in r 
when Of > 1. They have used the pointwise inequality (11.111) to establish this monotonicity formula. For the 
case of m > 2, we show that Ia{x) is monotone nondecreasing in r when 

2s$'(s) 


(1.23) 


a> a 


■- inf 

s>0 


$(s) 


We call this a weak monotonicity formula since the constant a* must be greater than one, due to some 
general assumptions on $. To clarify this, define an auxiliary function h(s) := — 2$'(s)s -I-a$(s) in the light 
of (11.231) . Note that from assumptions on $, i.e. 2s$"(s) -I- $'(s) > 0 when s > 0 and h{0) = Q;<i>(0) = 0 
one can see that h'{s) = — [2s<I>"(s) -|- $'(s)] -I- (a — l)<l>'(s) is negative when a < 1. For certain functions 
$, one can get the Laplacian, the p-Laplacian and the prescribed mean curvature operators and then the 
function Ia{r) is monotone in r when a > a* = 2, a > a* = p and a > a* = 2, respectively, see Corollary 
12.21 On the other hand, for both cases of scalar equations and system of equations, i.e. m > 1, it is shown 
in Theorem 12.41 that the following upper bound on the energy holds. 


(1.24) 


Ibr 


- 2H{u) + 2H{a) 


dx < CR'' 


where limx^^oo'Uiix' ,Xn) = at for all x = (x',Xn) S R" and a = This implies that for the case of 

system of equations, m > 2, the strong monotonicity formula, that is when a > 1, should hold just like in 
the case of scalar equations for m = 1. This remains as an open problem. Note also that conditions (A) and 
(B) are not necessary for our monotonicity formula when m > 2. 

We apply monotonicity formulae to establish Liouville theorems for solutions of (11.151) with a finite energy. 
We refer interested readers to Alikakos in and to Alikakos and Fusco in Q, to Caffarelli, Garofalo and 
Segala in 171 and to Farina in [s^, regarding Liouville theorems for various equations and systems with 
a finite energy. Note that the above monotonicity formulae, in both weak and strong forms, are related to 
the ones given for harmonic maps by Schoen and Uhlenbeck in and for minimal surfaces by Simon in 
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[5l|. by Ecker in [2^ and by Schoen in and for elliptic equations by Caffarelli and Lin in [18| , by Modica 
in [43l| and references therein. 

Regarding the scalar equation case, in this context, monotonicity of a solution u is straightforward to 
dehne and it refers to solutions that are monotone in one direction, e.g. when dx„Ui does not change sign, 
see iiilB , [ 2 ^, [ 23 , Isol - fs^ . [ 3 ^ . Issl 13,1131 and references therein. However, the notion of monotonicity of 
solutions for the case of system of equations, that is when to > 2, seems to be slightly more sophisticated. 
Ghoussoub and the author in [s^l introduced the following concept of monotonicity for the case of system of 
equations. Note that the sign of partial derivatives of the nonlinearity H could potentially have an impact 
on the monotonicity of solutions. This motivates us to call this notion as iL-monotonicity. 


Definition 1.2. A solution u = of is said to be H-monotone if the following holds, 

(i) For every 1 < i < m, each Ui is strictly monotone in the Xn-variable (i.e., dx„Ui ^ 0). 

(ii) For all i < j, we have 

(1.25) djHi{u)dx„Ui{x)dx„Uj{x) > 0 for all x G R”. 

See fsJ . [3/ for more details. 


Note also that in the assumption (ii) of the il-monotonicity each of djF[i(u), dx„Ui{x) and dx„Uj{x) 
has a fixed sign and the multiplication must be positive. This implies a combinatorial assumption on the 
system (11.151) . We refer to systems that admit such an assumption as orientable systems. For an example, 
consider to = 2 then for cross type solutions, i.e. dx„ui > 0 and dx„U 2 < 0, we are required to set 
diF[ 2 {u),d 2 Hi{u) < 0. If we set Hi{u) = H 2 {u) = — hu'i u'^ then this gives a two component system of 
equations that arrises in Bose-Einstein condensates, see ll| and references therein. 

The next definition is the notion of stable solutions for the case of system of equations. 


Definition 1.3. A solution u = (wfe)^]^ of U.15\) is called stable when there exists a sequence of functions 
4> = i4>k)'k=i such that each (jji does not change sign and djFli(u)(j)j(j)i > 0 for all i, j = 1, • • • , to. In addition, 
(j) satisfies the following 


(1.26) 


— div(.4(Vui)V(/)i) = ^ djHi{u)4>j in 


1=1 


where for any rj G M" the matrix A{r]) is defined by A{r]) := {a'i,jiv)) 7 j=i *11 U.KXl . 

Let us mention that the notion of stability can be given for weak solutions as 

/ m 

A{yu,)Vfi^ • VCi = ^ 

j=i 

where ( = (Ci)S=i is a sequence of test functions. Accordingly one can see that all results provided in the 
present paper are valid for weak solution as well. For the sake of simplicity in notation, we present our 
results for classical solutions. We refer to 21] and references therein for the use of stability for nonlinear 
elliptic eigenvalue problems. 

Since (|1.26|) is a linearization of (I1.15L one can see that every il-monotone solution is a stable solution 
via differentiating (11.1511 with respect to Xn- The notion of stability as well as the monotonicity formula for 
/^(r) when a > 1 and the pointwise inequality (11.1211 . provided by Modica, play key role in settling the De 
Giorgi’s conjecture (1978), see [ 2 ^. The conjecture states that bounded monotone solutions of Allen-Cahn 
equation are one-dimensional solutions at least up to eight dimensions. There is an affirmative answer to this 
conjecture for almost all dimensions. Moreprecisely, for two dimensions Ghoussoub and Gui in 37| and for 
three dimensions Ambrosio and Cabre in Q and with Alberti in Q gave a proof to this conjecture not only 
for Allen-Cahn equation but also for any equation of the form —Au = f(u) where / is a general nonlinearity 
that is locally Lipschitz. For dimensions 4 < n < 8 there are various partial results under certain extra 
(natural) assumptions on solutions by Ghoussoub and Gui in 38|, by Savin in and references therein. 


Note that there is an example by del Pino, Kowalczyk and Wei in [2^ showing that eight dimensions is the 
critical dimension. In two dimensions regarding the De Giorgi’s conjecture, we refer to Farina, Sciunzi and 
Valdinoci in [31 for a geometrical approach and to Modica and Mortola in [31 for some partial results under 
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the additional assumption that the level sets of solutions are the graphs of an equi-Lipschitzian family of 
functions. 

In 17|, I 22 I l33l| , authors considered quasilinear scalar equations of the form of (11.151) when m = 1 and 
provided one-dimensional symmetry and De Giorgi type results. Note that Ghoussoub and the author in 
[351 provided De Giorgi type results for elliptic systems of the form —Aui = diH{u) in lower dimensions for 
a general nonlinearity H. 

In this paper, we first provide a geometric Poincare inequality and a linear Liouville theorem for stable 
and iJ-monotone solutions of the quasilinear system (11.151) . Then we apply these to conclude De Giorgi type 
results for iJ-monotone and stable solutions in two and three dimensions when the system is symmetric. 
For coupled systems, that is when not all djHi vanish for 1 < i < j < m, it is natural to expect that there 
should be a relation between two arbitrary components Ui and Uj. In this regard, we show that gradients of 

all components of solutions are parallel and the angle in between and Vuj is precisely arccos ^ ^ 


when djHi ^ 0. This is a consequence of the geometric Poincare inequality, see Theorem 13.II 

The main focus of the present paper is the study of qualitative properties of solutions of system (jl.l5l) 
with a general nonlinearity. In this paper, we prove a Liouville theorem for bounded stable solutions of (jl.l5l) 
in dimensions n < 4 for a general nonlinearity H = {Hi)^^ whenever each Hi is nonnegative. To do so, we 
suppose that 4) satisfies either condition (A) or (B). Note that for the case of semilinear equations similar 
results are given by Dupaigne and Farina in and for the case of semilinear systems by Ghoussoub and 
the author in . In addition, we give a classification of radial stable solutions of symmetric system (11.151) 
when $(s) = |s 2 for all m > 1. More precisely, we show that there exists a positive constant Cn,m,p such 
that for any r, the following pointwise lower bound holds. 


(1.28) 


E 

2=1 


Uiir)\ > C, 


n,m,p 



ifn^ ^ 
ifn = ^ 

p-i 


-P, 

-P- 


This in particular implies that bounded radial stable solutions must be constant in dimensions 1 < n < 
-|- p. The notion of symmetric systems seems to be essential to study (ll.lSp with a general nonlinearity. 
Note also that the critical dimension n = -|- p for radial solutions is much higher than the dimension 

n = 4 derived for not necessarily radial solutions. Let us mention that for the case of semilinear equations, 
that is $(s) = s and to = 1, it is proved by Gabre-Capella 14, 15| and Villegas that any bounded radial 
stable solution of (11.151) has to be constant provided 1 < n < 10 when H G (^^(K) is a general nonlinearity. 
In addition, for the case of scalar equation and when $(s) = a counterpart of the above Liouvillle 

theorem is provided in 16, [T^. 

Here is how this paper is organized. Shortly after, in Section [5] we provide a Hamiltonian identity for 
solutions of system (11.151) . We also prove monotonicity formulae and we apply it to establish a Liouville 
theorem for solutions with finite energy. A few open problems are provided in this section as well. Section 
is devoted to some estimates needed to prove De Giorgi type results and Liouville theorems in next sections. 
We start the section with a stability inequality and then we apply this inequality to establish a geometric 
Poincare inequality. In Section SI we establish De Giorgi type results for iL-monotone and stable solutions 
of symmetric system (11.151) . In addition, we apply the geometric Poincare inequality, provided in Section 
m to find a relation between gradients of all components of solutions of ()1.15l) . Finally in Section!^ we 
prove Liouville theorems for stable solutions of (I1.15|) with a general nonlinearity, with on case requiring the 
solutions to be also radial. The concept of symmetric systems seems to be crucial to prove such an optimal 
Liouville theorem for radial solutions. 


2. Hamiltonian identities and monotonicity formulae 

We start this section by the following Hamiltonian identity. 

Theorem 2.1. Suppose that u = (uz)lLi is a solution of and let x = (x',Xn) G 1 

exists a constant C such that the following Hamiltonian identity holds for every Xn G . 


( 2 . 1 ) 


E(ivu,n-$'(|vu,nia,„u, 


X R. Then there 


- H{u) dx' = C, 
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when the above integral is finite for at least one value of Xn and in addition the integral in \2.1(]\] below tends 
to zero as R goes to infinity along a seguence. 

Proof. Suppose that x = {x',Xn) S K” and assume that -Bfi(O) is a ball of radius R in Define 

r : R —>■ K as 


( 2 . 2 ) := 

Differentiating T with respect to Xn we get 


/ E 

IbrIo) Vfci L 


- H{u) ) dx'. 


f 1 

-Hi{u)dx„Ut}dx' 

m « 

(2.3) =: E / {ri(cc)+r"(a:)+r3(x)}dx'. 

j_i JBr{o) 

In what follows we simplify the above three terms, appeared in the right-hand side of (12.3|) . Note that 

(2.4) [$(|Vu,n] = 2$'(|Vu,nVu, • 

Therefore, 

(2.5) ri(:c) = $'(|Vu,nvu, • 

Similarly, 

(2.6) r2(x) = -2$"(|vu,nvu, 

-2^' {\yui\'^)d^^u^dY^u^. 

We now apply (11.1511 to simplify T^, 

r^(a;) = div($'(|V?.tip)Vui)93,„Ui = div^,-($'(|Vuip)Vx'Ui) 

+2$" (|Vu,p) Vu, ■ 

(2.7) 

Adding (12.5|) . (12.61) and (12.71) we get 

(2.8) ri(a;)+r2(x)+r3(x) = div,, ($'(|vu,nv,m,)a,„u, 

-f$' (|Vuip) d^^UidY^Ui. 

Substituting (12.81) in (12.311 and applying the divergence theorem we obtain 


(2.9) Tfiixn) = V/ 

Jobr^o) 

Suppose that the integral in (12.11) is finite when Xn = 0. Then, 

( 2 . 10 ) TR{xn)-rR{0) = Y / 

Jo JdBRiO) 

Taking the limit of the above when i? —?> oo finishes the proof. 

□ 


When $ is the identity function, the system of equations (11.1511 is a semilinear system of the following 
form 

—Aui = Hi{u) in R”. 

Therefore, Theorem 12.II implies that the following Hamiltonian identity holds. 



dx' = C. 
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( 2 . 11 ) 










Note that (12.111) is given by Gui in . Here we have Hamiltonian identities for the mean curvature system 
as well as the p-Laplacian system. 

Corollary 2.1. Suppose that assumptions of Theorem \2. 1\ hold. 

(i) Let d)(s) = 2{^/n^ — 1). Then U.15\) reads 

Vui 


— div 


yi + ivM,|2 

For any Xn G K., this Hamiltonian identity holds, 


= Hi(u) ir. 


( 2 . 12 ) 


/B"-i 


^i=i 


1 + iVx'Uip — y^l + |Vui| 

\/l + |Vu*P 


- H{u) dx' = C. 


(ii) Let $(s) = |s 2 . Then i f 1.151) reads 

— div (^\'\/ui\^~^'Vui) = Hi{u) ir, 
For any a;„ G M, this Hamiltonian identity holds, 


(2.13) 


/R"-i 




P-2 


P p 


- H{u) dx' = C. 


where C is a constant. 


Note that for a specific nonlinearity of the form Hi{u) = m [l — J2rLi original system of equations 

(fTTKl) is 

— div($'(|VMip)Vui) = Mi ^1 — ^in R". 

This is a quasilinear Ginzburg-Landau system where u = {ui)'fLi for Ui : K" —>■ K. For the semilinear case 
and m = 2, see [l3)ll3- Since H = — YllLi ^f)^ is an antiderivative of H, the following Hamiltonian 

identity holds, as long as conditions of Theorem 12.II are satisfied, for any G M 

/ ™ \ 2 ' 


(2.14) 


E 


E(|VM,|2)-$'(|VM,|2)|a,„M,|2 


1 


ih-E' 


dx' = C. 


For the rest of this section we study monotonicity formulae for solutions of system (I1.15|) . Gonsider the 
following function (r) for any r > 0 


(2.15) 


!a(r) := 


^$(|Vm,|2)-2H(m) 


i=l 


For the case of scalar equation, that is when m = 1, the following monotonicity formula holds. Note that 
this is a direct consequence of Theorem A. 

Theorem B. 0/ Suppose that m = 1 and u is a solution of djf). In addition, suppose that assumptions 
of Theorem A hold. Then the functional Ia{r) when a>\ is a monotone nondecreasing function of r. 

For the case of system of equations, that is when m > 2, we provide a weaker version of the monotonicity 
formula provided in Theorem B, under certain lower bounds on a depending on <i>. 

Theorem 2.2. Let u = (ui)(T]^ be a solution of ill.151) when m >2 and H{u) < 0. Suppose that there exists 
a constant a such that 


(2.16) 


^ , . J2s<l>'(s)) 

a>a := mf 1 , > . 

s>0 ( <i)(s) 


Then, the functional lair) is a monotone nondecreasing function of r. In particular. 


(2.17) 




[ f;$'(|VM.n(a.My--JET/ Hiu). 

J dBr i —\ J Br - 

































Unlike Theorem B, conditions (A) and (B) do not appear in assumptions of the above theorem for the 
case of system of equations. This implies that Theorem 12.21 is valid for a larger class of nonlinearities $ 
compared to Theorem B. However, as mentioned in Section [1] the constant a* must be greater than one, 
due to assumptions on $. For the sake of convenience of readers we clarify this here as well. Consider the 
auxiliary function h{s) := —2<i>'(s)s + a$(s) regarding terms appeared in (12.161) . Note that from assumptions 
on $, i.e. 2s<i>"(s) + d>'(s) > 0 when s > 0 and h(0) = a$(0) = 0 one can see that h'{s) = — [2s<i>"(s) + 
$'(s)] + (a — l)$'(s) is negative when a < 1. We now compute a*, provided in (I2.16L for various choices of 
$. 


Corollary 2.2. Suppose that m>2 and H <0. Consider the following particular functions $. 

(i) Let <f>(s) = s and u = (ui)™ i he a solution of the semilinear system of equations 

(2.18) -Au^ = Hi{u) in R”. 

Then for all a > a* = 2, the function lair) is monotone nondecreasing in r. 

(ii) Let <l>(s) = 2(\/l + s — 1) and u = (ui)™ i be a solution of the mean curvature system of equations 


(2.19) 


— div 


Vui 


\/l + |Vu*P 


= Hi{u) ir, 


Then for all a> a* = 2, the function lair) is monotone nondecreasing in r. 


(iii) Let d)(s) = and u = iui)'^^ be a solution of p-Laplacian system of equations 

(2.20) -div(|VM,|P-2VM,) =Ff,(u) mR". 

Then for all a > a* = p, the function lair) is monotone nondecreasing in r. 

To provide a proof for Theorem l2.2l we present a few technical estimates. We follow a classical argument 
regarding Pohozaev and Rellich type identities ^ to prove the following identity on a ball of radius r. 

Lemma 2.1. Suppose that u = (ui) is a solution of il.l5\) then 

-n V$(|Vu,|2) = 2r V $'(|Vu,|^)(9,u,)'- ^ 

Job, Job, 

n 77 % p p 

-2 ^$'(|Vu*n|Vu*p - 2n / Hiu) + 2r Hiu). 

B-r 2_ Br- 

Proof. Multiply the i*^ equation of (11.151) with x - S/ui and then apply the divergence theorem to get 

[ $(ivu,n = 2 / $'(|vu,n(a,u,)2 - ^/■ $'(|vu,nivu,p 

Job, Job, Jb, 

[ (x • Vui) div ($'(|VMip)Vui) + - / ^dVuip), 

^ Jb, ^ Jb, 

for each i = 1, • • • , m. Doing some straightforward computations as well as applying (11.151) gives the desired 
result. 

□ 


We are now ready to provide a proof for Theorem 12.21 
Proof of Theorem \2.2[ Differentiating /(r), given by (12.151) . with respect to r gives 


I'air) 


.,71 —a + 1 


p 771 

= ia-n) V$(|Vud")-2i?(u) 

p 771 

+r V$(|Vu,n-21?(u). 

JdB ,. 1 
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Substituting the value of —n S'S it is provided in Lemma I^TTl one can show that 

P m 

= 2r 

JdBr 

c r ^ 

+ / ^ (-2$'(|Vud")|Vu.p + a$(|Vud")) - 2aiL(u) 




i=l 


The rest of the proof is straightforward. 


□ 


As an immediate consequence of Theorem 12.21 we have the following Liouville theorem for solutions of 
(11.151) with a finite energy. 

Theorem 2.3. Suppose that assumptions of Theorem \2.2\ hold. Assume also that u = (rtz)™i has a finite 
energy that is 


( 2 . 21 ) 




i=l 


dx < oo. 


Then each ut must he constant in dimensions n > a for i = 1, ■■ ■ ,m. 

Proof. First suppose that n > a. From Theorem 12.21 for any R> r we have 

1 


0</„(r) = 


< 


< 


i?" 


Ibr 


i?" 


m 

m 

^$(|Vu.n-2iL(u) 
.2 = 1 


Sending i? —^ oo, we get the desired result. Now suppose that n = a. Again from Theorem 12.21 we have 

(2.22) rl'^ir) > -2a [ iL(u). 

JBr 

From this for any r > f, where r is fixed, we have 

(2.23) Ia{r) > Ia{r) —‘2a\n j H{u). 

Note that (12.211) implies that limr_>.oo Ia{r) < oo. From this and (12.231) we conclude that H = 0. The fact 
that each component m is harmonic together with (12.211) completes the proof. 

□ 

Another consequence of the monotonicity formula, given in (12.171) . is the following lower bound on the 
energy. 


Corollary 2.3. Suppose that assumptions of Theorem \2.A hold. Then, the following lower bound holds for 
the energy functional 


(2.24) 


'Br 


5]$(|Vu.n-2iL(u) 


. 2=1 


dx > CR°‘-^ for all R>1, 


where C = 1(1) is independent from R. 

As the next theorem, we prove an upper bound on the energy function. 
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Theorem 2.4. Suppose that u = (ui)^^ is a bounded H-monotone solution of such that for each 

i = 1,m, 

(2.25) lim Ui{x',Xn) = cii, Va: = (cc',a;„) G M" 


for some constants ai. Then 
(2.26) Jr{u) := 


ibr 


^(iVitip) - 2H{u) + 2H{a) 


dx < CR” 


where a = (az)i(Li C' are independent from R. 

Proof. Define the sequence of shift functions m* = (m*)™! where ul{x) := Ui{x',Xn + t) for t G M and 
X = (x',Xn) G M". Note that u* = {u\)'ffi satisfies 

(2.27) -div($'(|VM*nVu‘) =i7*(u‘) in R’". 

The fact that u\ is convergent to Oi pointwise, it is straightforward to see that 


(2.28) 


'Br 


{H{u*) — H{a))dx ^ 0 when t —?> oo. 


On the other hand, multiply both sides of (12.271) with u- — at and integrate by parts to get 

- / <i>'(|V«r)|V«‘p + / $'(|v«r)a.«‘(u‘ -«,) = -/ - «.), 

Jbr Jobr Jbr 


I Br 

Sending t —?> oo implies that 


[ $'(iv«‘nivu‘p^o. 

Jbr 


Note that due to the assumption 2s$"(s) + $(s) > 0 when s > 0 and $(0) = 0 we have 0 < $(s) < 2$'(s)s 
for any s > 0. This implies that 


(2.29) 


0< / $(|Vu‘n<2 [ $'(|Vu*n|Vu*|^ ^ 0 as t ^ 0. 

JBr JBr 


From this and (12.281) we get 
(2.30) 


lim Jr{u*) = 0. 


We now use Jr{u*) to construct an upper bound on Jr{u). Note that differentiating the energy functional 
with respect to t, one gets 

m « 

(2.31) dtJRiu*) = V / [2<l>'(|Vu*nVu‘ • Vidtul) - 2iJ,(u‘)atu*] . 

i=l -^Br 

Multiplying the system of equations (12.271) with dtvf and performing integration by parts we obtain 

(2.32) f m{u^)dtu\ = f $'(|Vu‘^Vu^V(a^un 

JBr Jbr 


IOBr 




for each i = 1, - ■ ■ , m. From (12.321) and (12.311) we obtain 

m p 

(2.33) dtER{u*)=y^ ^'{\Vu*f)dRuldtul 

i=l JoBr 

Note that there exist a constant M such that —M < $'(|Vu-p)9,yU* < M and dtu* > 0 > dtu* for i G / and 
j G J ■ Therefore, 


(2.34) 


dtERiu*) > M dtul j dS. 


j&J i€l 
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Therefore, 


Jr{u) = Jr{u*) - f dtJR{u‘')ds 
Jo 


^ Jr{v^^ T -Tf 


daUl — daU’^j dSds 

0 


(2.35) 


= Jr(u^) + M [ [ y^(M- - Ui) + y^(Mj - m‘) I dS. 

V.e/ / 


From the definiton of iJ-monotonicity and the sets of I, J we have ut < u\ and m*- < Uj for all i G I, j € J 
and t G M"*". Therefore, 

(2.36) Jr{u) < Jr{u*) + M j I y^(Mi — Ui) + y^(uj — u*) I dS for all t G K'*', 

JdBR \ T / 


i€l 


j&J 


where M := max™^ {||d>'(|Vuip)|Vui|||Loo(]an)}. The upper bound (12.361) implies 

Er{u) < Er{u*) + C\dBR\ for all t G K'*'. 

Sending t —?> oo and using (12.301) . finally we obtain that 

jR{u)<C\dBR\<CR^-\ 

This provides the desired result. 


□ 


Before we finish this section, let us mention a couple of open problems for the system of equations (11.151) . 

Open Problem 1. Under what assumptions on H = (Hi)^-^ and solutions, one can provide a eounterpart 
of the pointwise inequalities provided by Modica in fISl and Cajfarelli et al. in fT^I for solutions of \1.15\) 
when m > 2? 

Open Problem 2. In the light of Theorem \2.4\ Corollarv \2.S\ and Theorem B, one might expect that lair) 
should be a nondecreasing function of r when a > a* = 1 for the case of systems that is when m >2. 

3. Geometric Poincare and stability inequalities for systems 
Note that the matrix A{r]) := (aij(^))Lj=i where aijig) is defined by 

(3.1) a.,M■■=^•^”i\d\^)v^rlJ+^'i\rl\^)S^J■ 

is symmetric and positive definite for every rj G K". This is because for any f G M", 


(3.2) 


Aiv)c.c = E 2c*o‘i>"(bnw,+‘&'(bnGo<5.., 

= 2$"(|,yn|C-r;p + $'(|^|2)|^|2, 


Note that when ^"(rj) is positive clearly Air])(4.( is positive since is positive and when $"( 17 ) is negative 
applying Young’s inequality together with 2$"(s)s + $'(s) > 0 when s > 0 implies that Airi)f.() is positive. 
We are now ready to prove the stability inequality for solutions of (11.151) . Note that such an inequality for 
the case of semilinear systems is given in [^, 1 ^, . 

Lemma 3.1. Let u = (ui)™ denote a stable solution of il.l5\) . Then 

(3.3) y] / JdjH,iu)d,H,iu)QCj <Y. -4(Vw,)VG • V^, 

for any C = '^here (i € for 1 < i < m. 


12 















Proof. Since m is a stable solutions, there exists a sequence (j) = that satisfies (ll.26p . Consider a test 

function ( = where Q G fl with compact support and multiply both sides of (11.261) 

with Integrating by parts we get 


(3.4) 




A(Vui)'V(pi 




=.Ji. 


The fact that v4(Vui) is positive definite we get 


0 < A{yUi){(j)iS/C,i - • {(fiS/Qi - CiS/fn) 

= (flA{Vui)VC,i ■ VCi + C,];A{Vui)V4>i ■ - Q(j)iA{\7Ui)V(j)^ ■ VC*. 


Applying this to (laa for each i we obtain 


(3.5) 

For the left-hand side we have, 


< [ A{Vu^)V(i-VCi. 


I 




jri . ^2 ^ p p2 p 

Y, +Y Wu)cj,,f + Y 


i<j 


i>j 


i—1 ' 

^2 


> 


Tfl n ^'2 '7* p p2 p 

Y d,H,{u^,f + Y d.H,{u)cf,f+Y 

m „ / .2 \ ™ /• 

^ )+Y 

i<j 7*" \ JR" 

m « _ m « 

2Y Jd,m{u)mxu)QQ+ Y / 

7Rn Jr- 

Y f \/djH,{u)d^HXu)CiCj. 

*, 7=1 7k" 


This finishes the proof. 


□ 


We now apply the stability inequality to provide a geometric Poincare inequality of the following from. 
For the case of scalar equations that is when m = 1 this inequality was driven by Sternberg-Zumbrun in 
( 5 ^ and it was applied in this context by Farina-Sciunzi-Valdinoci 331 and references therein to provide 
De Giorgi type results. Note also that Cabre applied this type inequality to prove regularity of extremal 


solutions of nonlinear eigenvalue problems in 


3 app 

0. 


For the case of system of equations that is when m > 1 


this inequality was first proved by Ghoussoub and the author in 
conclude De Giorgi type results for system of equations. Let us mention that interested readers can find 
similar geometric Poincare inequalities in these references as well 


25, 


and they applied the inequality to 
on tf 

Mim, 


13 













Theorem 3.1. Assume that m,n > 1 and u = (ui)™i is a stable solution of Then, for any 

V — iVk)'k=i £ C'g (R"), the following inequality holds; 


(3.6) 


/ JdjHi{u)diHj{u)\Wui\\Wuj\riir]j - djHi{u)Wui ■ WujT]^ 

/tC)71 ’ 




+ / <i>'(|Vu,|2)|Vu,|2/c2r,f 

m p 

+ E / [2$"(|Vu,n|Vu?| + $'(|Vu,n] |Vt.|Vw,||V 

Z—1 

m p 

< V / 2|Vu,|2$"(|VM,n|VM, •Vr,,p + $'(|Vu,n|Vu,nVry,p, 

JR" 


i=l 


where Vt^ stands for the tangential gradient along a given level set of Ui and Kf for the sum of squares of 
principal curvatures of such a level set. 

Proof. Suppose that 77 = ( 771 , 77 ^) for 77 ^ € is a test function. Test the stability inequality (15.81) 

with fi = |V77,i|77i to get 


(3.7) 



dtH,{u)\\7ui\'^Pi 


< 




|Vuip.4(VMj)V77i • \7ri^ 


m p 

+ E / »7*^-4(V'Ui)V|V77j| • V|V77j| 

i=l"" 


- 5 I: 

2=1 ‘ 


.4(Vu,)V|V77,|fo V 772 . 


Straightforward calculations show that for each k, 

(3.8) d,, ($'(|V77,nV7i,) = .4(V7X,)Va,,77,. 

Applying this and differentiating the 7*^ equation of (I1.15|) with respect to Xk for each i = 1,2, ...,m we get 

m 

(3.9) - div {A{Vui)\7dxkUi) = Y.diH.{u)dx,u.. 

i=i 

Multiplying the above with gfdkUi, integrating by parts and taking sum on the indices i, k we get 


(3.10) 


m p 

Y, / a.i7.(u)|Vu.p77f 

i=i “'R" 

P m n 

- E / 9jHi{u)Wui ■ Vujpi + E E 

i=i k=i 

P m n 

- E / 9jHi{u)Wui ■ Vujpj + E E 

m „ 

+ oE / - 4 (V 7 x,)V|Vu,|foV 772 . 


2 = 1 /c=l ' 


A{yu,)v{dx,uf) ■ v{dx,v!) 


77,M(Vui)V(a^^'Uj) • V{dx,,Ui) 
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Equating (13.101) and (13.71) we get the following since the term i Y^=i /r>* ^(Vui)V|VMip • S/rj^ cancels out, 


(3.11) 

(3.12) 

(3.13) 

(3.14) 


/ J djHi{u)diHj{u)\'S/ui\\Vuj\'qi'qj - djHi(u)Vui ■S/Uj'ql 

l 

m « 

+ / y4(Vui)V(5a;;,'Ui) • 

m „ 

-V / ^(Vu,)V|Vu,| • VlVu^lry^ 

i=i -'R" 

m ^ 

< ^ / |Vu*p^(Viii)Vry, ■ V/7i. 


For the rest of the proof, we simplify two terms (13.121) and (13.131) in the left-hand side of the above inequality. 
From the definition of there matrix A, in the light of (13.21) . we get 


(3.15) y]^(Vu,)V(5,,u,)-V(5,,zi,) = 2$"(|Vu,|") |Vu, • 

k^l k^l 

n 

+4>'(|vu,ny]|va,,u.p. 

k^l 

Straightforward calculations show that 

n 

(3.16) y] |Vu, • = - |V|Vu,p|" = |Vu,p |V|Vu,||" . 

A.-1 

From fl3.15p and (|3.16p we obtain the following form for the term in (|3.12p 

m n p 

i=i fe=i -'R" 

m p 

(3.17) = Y 2|Vu,p$"(|Vu,n|V|Vu,||'r?f 

i=i -^R” 

m p n 

+ Y <i>'(|vu,ny]|va,,u.|V- 

JR”- 


k=l 


Similarly, from the definition of the matrix A, i.e. using (13.21) . we get 

^(Vu^VlVuil • VlVuil = 2$"(|Vu,n |Vu, • V|Vui||^ 

(3.18) +$'(|Vu,n|V|Vu,||". 

This implies that the term in (13.131) is of the from 

m p 

Y ^(Vu,)V|Vu,|-V|Vu,|7?f 

i=i “'R'‘ 

m p 

(3.19) = E / 2$"(|Vu,n|Vu, • V|Vu,||"7?2 

i=i "'R" 

m p 

+Y $'(ivu.nivivu.i|2 7?2. 
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The difference of (13.171) and (13.151) . as appeared in (13.131) and (13.121) . is 


(3.20) 

(3.21) 

(3.22) 

(3.23) 


m n ^ 

^ ^ / A{Vu.)V{d,,Ui) ■ 

i=l k=l 
m ^ 

-V / .A(Vu,)V|Vn,| • V|Vn,|??2 


= 2j2 / 




|V|Vmi||^ - —-—|VMi • V|Vui 


m « 

+ ^ / $'(|Vw.|2) 
,^iJn 


jVui 

n 

- \V\Vu,\ 


$"(ivw,n 


lk=l 




where fl = {|VMi| ^ 0} (h R.". We now simplify (13.221) and (I3.23|) via applying the tangential gradient and 
curvatures. Suppose that \Vui \ ^ 0 at a point x G M”, then 


(3.24) 


|Vt|Vm.||2 = |v|Vm.|| - 


1 


\VUi 


|Vmi-V|Vmi|| , 


where Vr denotes the orthogonal projection of the gradient along this level set. In addition, according to 
formula (2.1) given in the following geometric identity between the tangential gradients and curvatures 
holds, 


(3.25) 


^ - |V|Vn,|p = |Vn,|2/C2 + |Vt|Vu 




fc=i 


for /Cf := where ki is the principal curvatures of the level set of Ui at x. Substituting (13.241) and 

(I3.25|) in (13.221) and (13.231) we get 


(3.26) 


m n p 

i=l k=l 

m p 

-Y. A{Wu,)W\Wu,\-W\Wu,\r]f 

.1 JR-^ 


= E 


2|VM,|"$"(|VM,|")|Vr|Vu,|p772 


_i {IVuiiT^ojnM.”' 


+ ^ / $'(|Vw.n [|Vu,p + |VT|Vn,||^] r,l 

i—l -^{l 

Finally, substitution of (I3.26|) in p.l3l) and (13.121) completes the proof. 


□ 


4. De GiORGI type results for SYIMMETRIC SYSTEMS 

In this section, we provide One dimensional symmetry results for stable and iJ-monotone solutions of 
symmetric system (11.151) in lower dimensions with a general nonlinearity. At first let us fix a few notations. 
Throughout this section we suppose that C, = (Ci)i2.i is a sequence of test functions where Q G C'^(R"’) (~) [0,1] 
where C* = 1 in Bi and C* = 0 in R"- \ B 2 . Set rQ{x) := 0(|) and /jT = (^X)™ ^ where i{ri(a:) := VCj(f) 
for any i? > 1. Note that ||fl^i(a;)|| 2 ,oo(R,^) < C where C is independent from R and S/RQi{x) = R~^RTi{x). 

To set up a Liouville theorem for the quotient of partial derivatives of solutions of (11.151) . we first state 
the following technical lemma. 
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Lemma 4.1. Suppose thatu = is aH-monotone solution of il.l5\) . Set(j)i := dx^Ui andifji := ’S/ui-rj 

where p = (p',0) G x {0} and define Ui := Then the sequence of functions a = (cTi)™ satisfies 

m 

(4.1) div [(j)iA{Wui)Wai\ +'^djHi{u)(l)i(j>j{aj — ai)ai = 0 in R”. 

1=1 

Proof. Since the proof is straightforward we omit it here. □ 

The fact that a = (ci)™ ^ satisfies 631 motivates us to provide a Liouville theorem for system dm. 
Applying Caccioppoli type arguments we establish the following Liouvlle theorem for a slightly more general 
setting than m- Let us mention that for the case of scalar semilinear equation, m = 1 and 4>(s) = s, this 
type of Liouville theorem was noted by Berestycki, Caffarelli and Nirenberg in and used by Ghoussoub 
and Gui in 371 and later by Ambrosio and Cabre in @ to prove the De Giorgi conjecture in dimensions two 


and three. Also, Ghoussoub and Gui in [38| used a slightly stronger version to show that the De Giorgi’s 


conjecture is true in dimensions four and five for a special class of solutions that satisfy an antisymmetry 
condition. We also refer interested readers to Q by Barlow, Bass and Gui and to @ by Barlow for some 
probability based arguments regarding this Liouvlle theorem. 

For the case of scalar quasilinear equation, m = 1 and a general $, this Liouville theorem is provided by 
Farina, Sciunzi and Valdinoci in 331 and by Danielli and Garofalo in (2^. For the case of semilinear system 
of equations, m > 1 and $(s) = s, we refer to [S^ by Ghoussoub and the author. 


Proposition 4.1. Assume that for each i = 1, 
where (fl > d and Oi G Let 


(4.2) 


E 


B2r\Br 


,m functions |Vu^| and (pi are locally hounded in 


where the constant C is independent from R > 1. Let a = (cr^)™ satisfy 

m 

(4.3) (Ti div [0^A(Vui)Vcri] + ^ h^j{x)f{aj - ai)a^ >0 in R", 

i=i 

where 0 < hij G Lj^q^(R"), hij = hji and f G L)g^(R) is an odd function such that f(s) > 0 for s G R’*'. 
Then, each function Ui is constant for all i = 1, ..., m. 

Proof. The proof is strongly motived by the ideas and methods used in ii [ssl - l^ . Note that 

m 

hij{x)aif{aj - af) = hijUifiaj - cri) + hijaif{<7j - at) 

i,j=l i<j i>j 

= E + E ~ since % = hji 

i<j i<j 

= — hiji^j — <^i)f{<^j — <^i) since / is odd. 

i<j 

< 0 since hij > 0 and s/(s) > 0 for any s G M. 

Multiply (14.31) with a text function PiC,i^ and perform integration by parts to get 

m « m « 

v/ 0iA(Vuj)V(CTji{C^) ■ VtTj < V / htj{x)aJ{aj - a^) <0. 




Therefore, 

(4.4) 


E / ■ Vcr* < -2^^ 


i=l ' 


— Id B2 r\Br 


■ VcTi. 


In the light of the Gauchy inequality with epsilon, one can see that for any e > 0, 

2\ai\H(i\A{\7ui)VH(i ■ VcTil < A(VMi)Vcri • Vcr* + -afA{Vui)\7 rQ ■ nfi. 

e 
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From this and dUl), for every R > 1, we get 


(4.5) V/ • Va, < • Va, 

1 m 

H—^ / <('i-4(Vui)V_RCj ■ V_rCj. 

^ i=l B2r\Br 

From this, (j4.2D and the definition of the test function nQ we conclude that the following integral is bounded, 


(4-6) E / (/)f.4(Vui)V(Ti • Vdi < oo. 

i=l "'k" 

Now, sending R ^ oo and e —>■ oo and applying (H3|) and (14.21) show that the integral in (14.61) vanishes 
for every i = I,-- - ,m. Finally, the fact that .4 is a positive definite matrix implies |Vcri| = 0 for each 
i = 1, - ■ ■ , TO. This completes the proof. 

□ 


We are now ready to present the following De Giorgi type results in two dimensions. 


Theorem 4.1. Suppose thatu = is a classical bounded stable solution of symmetric system S1.15\) in 

two dimensions. Assume also that |Vui| G (R^). Then each Ui is a one dimensional function, 

i.e. there exists u* : K —>■ R and a G such that Ui(x) = u*(a ■ x). In addition, the angle between "S/ui and 


Vuj is arccos 


V SiHi ) ■ 


Proof. We apply the geometric Poincare inequality given as Theorem 13.11 to provide a proof. Ideas and 
method applied in this proof are strongly motivated by the ones given for the case of the scalar equation 
by Berestycki, Caffarelli and Nirenberg in Ghoussoub and Gui in [s^, Farina, Sciunzi and Valdinoci in 


331 and references therein. In addition, for the case of system of equations we refer interested readers to 
35 ,( 3 ^ by Ghoussoub, Sire and the author. Note that from boundedness of |VMi| and G C'(R+) in two 


dimensions we have 
(4.7) 


IBr 


$'(|Vu,n|Vu,|2 < CR^ 


for any R > 1. This can be also proved by multiplying (11.151) by nCfui and integrating by parts. Straight¬ 
forward calculations show that for each i we have 


1 


I Br\B^ 




< 2 


2 / [ T-3$'(|Vu.n|Vu.|2drcix + ^ / $'(|Vu,n|Vu,| 

Jbr\b^J\x\ ^ JbrXb^ 

r T-^ [ $'(|Vu,n|Vu,pd2:dT + 4 / $'(|Vu,n|Vu,p, 

Jvr Jb., R Jbr 


where we have used the Fubini’s theorem. From this and Kn we get 

f $'(|VU,|2)|VU,|2 


(4.8) 


< Clogi?. 


JBr\b^ fol 

Now for each i set rji to be the following standard test function 


1 

2 ’ 


if \x\ < \/R, 

r,,(x):={ H VR < \x\ < R, 

0, if |a;| > R. 
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Note that Vrji{x) = — ^ 1^1 ^ Therefore, the right-hand side of the inequality given in 

theorem o is of the form 


< 


< 


< 


m „ 

V / 2\Vu^f<^>"{\Vu^f)\Wu^■Vr]^f+ <^>'{\Vu,\^)\Wu,\'^\Vf]^f 

i=i 


m « 

/ \\/ui\'^A{Vui)Vf]i-Vrji 

j_l JBr\B^ 


C 

log^R 


i=l ''Br\B^ 


C 

\og^R 


E 


= 1 J Br\B^ 


Ui\^A{VUi)x ■ X 
\x^\ 

$'(|Vn,|2)|Vn,|2 

It-|2 


c 

logi?’ 


where we have used (14.81) to conclude the last inequality. From this and the geometric inequality given by 
dSH) and the fact that the system is symmetric we get 


(4.9) 


V [ {\djH,iu)\\Vu,\\yuj\- d,H,{u)Vu, ■ Vuj) 




E / $'(|Vn.|2)|Vn.|2/Cf 

j_i J{\vui\^^}r\B 


E 


< 


i=i •^{|Vui|^o}nB^ 

C 


[2$"(|Vu,n|Vu2|-f $'(|Vu,n] |VtJVu,| 


logi?' 


Now sending i? —>■ oo and the fact that all of the terms in the left-hand side are nonnegative imply that each 
Ui is one dimensional function and \djHi{u)\\S/ui\\S/Uj\ = djHi{u)'S/ui ■ 'S/Uj. The latter implies that the 

angle between Vui and Vuj is precisely arccos when i ^ j. This completes the proof. □ 


Note that in the statement of Theorem EH $ does not have to satisfy conditions (A) or (B). However 
in the next theorem that is in regards to three dimensions one of conditions (A) or (B) is needed. For the 
case of semilinear systems in two dimensions we refer interested readers to [l| for the construction of two 
dimensional solutions in the absence of stability and iJ-monotonicity, 

Theorem 4.2. Suppose that u = (ui)™! is a classical bounded H-monotone solution of symmetric system 
m.l5\) in three dimensions. Let $ satisfy one of conditions (A) or (B). Assume also that |Vrti| £ L°°(]R^) D 
Then each Ui is a one dimensional function. In addition, the angle between Vui and Vuj is 
arccos • 

Proof. Methods and ideas applied here are strongly motived by the ones given by Ambrosio and Cabre in @ 
and Alberti, Ambrosio and Cabre in in the case of a single equation and by Ghoussoub and the author 
in for the case of systems. We first note that u being Ff-monotone means that u is a stable solution 
of (11.151) . Moreover, the function Vi(xi,X 2 ) := linix^j^oo Ui(xi, X 2 , X 3 ) is also a bounded stable solution for 
(jl.lhp in K^. Note also that since u is an iF-monotone solution, the system (11.1511 is then orientable. It 
follows from Theorem 14.11 that each Vi is one dimensional and consequently the energy of u = {vi)'fLi in a 
two-dimensional ball of radius R is bounded by a multiple of R which implies 

(4.10) limsupF;(u*) < Ci?2, 

t—¥oo 
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where u*{x') := u{x', Xn + i) for t G K and 

P m 

Er{u)= / V$(|V«,n-H(M) + c„, 

■IBr 

for Cu := maxi7(rt). Applying similar arguments as in the proof of Theorem l2.4[ we shall show that 

m p 

(4.11) Y. nWu,\^)<CR\ 

i=i Jbr 

where the constant C is independent from R. Note that shift function u* = (m*)™ i is also a bounded solution 
of (fTTKll with |Vm*| g i.e., 

(4.12) -div($'(|Vu‘nVM*) =i?i(w*) in R", 
and also 

(4.13) dtu* > 0 > dtu* for alH G / and j G J and in R". 

Since u- converges to Vi in C'/q^(R"') for alH = 1, • • • , m, we have 

lim E{u*) = E{v). 

t—¥oo 

Now, we claim that the following upper bound for the energy holds, for all t G R"*" 


(4.14) 


Er{u) < Enin*) + M f [ ^(m- - Ui) + y'(uj - u‘) | dS, 


where M := max™;^ {||$'(|Vuip)|Vui|||icx=(Rn)}. Indeed, by differentiating the energy functional along the 
path M*, one gets 

m p 

(4.15) dtER{u^) = Y 2$'(|Vu‘nVu‘-V(atu‘)-2i?,(w*)atM^ 

i=l JBr 

Now, multiply (14.121) with dtu^ and integrate by parts to get 


(4.16) 


'Br 


+ [ ^'{\\7ul\'^)d^uldtul = [ Hi{u^)dtu\. 

JOBr JBr 


for each i = 1, • • • , m. From (14.161) and (14.151) we obtain 

(4.17) dtER{u^) = 2Yl ^'{\^u\\^)d,u\dtu\. 

i JdBR 

Note that —M < $'(|Vu-p)9i/M* < M for each i and dtu\ > 0 > dtvd, for i G / and j G J ■ Therefore, 


(4.18) 

On the other hand. 


dtERiu*) > M 


/OBr 


YdtuYY.dtul]dS. 
\jeJ iei 


Er(u) = Er{u*) - [ dtER(u‘^)ds, 
JQ 


< Er{u*)+MI I I Y, ^sUj - I dSds 

jei jeJ 


0 JdBR 


(4.19) 


= ER{u*)+Mf [ W(m-- Mi) + - u‘) I dS”. 
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To finish the proof of the theorem just note that Ui < u\ and u* < Uj for all i G I, j G J and t G K+. 
Moreover, from (14.101) we have Eii{u^) < CE?. Therefore, (I4.19P yields E}i{u) < C\dB}i\ < CR^. 

This proves (14.111) . 

Now, set (j)i '■= dx^Ui and := Vui • rj where rj = {r]',0) G x {0} and define Ci := Lemma l4T] 
implies that tr satisfies 611. Note that ipfcrf = ipf < |Vuip. From this and the fact that one of conditions 
(A) or (B) holds there exists a constant M that is independent from R such that 


(4.20) 

(4.21) 


i=l ' 


B2 r\Br 


cl)faiA{\7ui)Rri ■ liTi 


771 p 

< ^ Ui) rT i ■ rT i 

i=l JB2R 

m p 

< 

i=l 


here we have used the fact that ||_Rrj(a:)||ioo(Rn), ||VMi||j;,oo(RTi) < C for some C that is independent from R. 
We now apply Proposition 14.II for /i^ = djE[i{u)4>i(j)j and the identity function / to conclude that each Ui 
is constant. This implies that each Ui is a one-dimensional function. Since u = (u^)™ is a one-dimensional 

stable solution of (I1.15L Theorem 14.11 implies that the angle between and Vitj is arccos 

□ 


5. Liouville theorems for symmetric systems 


For bounded stable solutions of (11.151) up to four dimensions we have the following Liouville theorem as 
long as each E[i(u) is nonnegative. 

Theorem 5.1. Suppose that u = (itz)™! is a classical bounded stable solution of symmetric system \1.15\) 
where iJi > 0 for each i. Let $ satisfy one of conditions (A) or (B). Assume also that |VMi| G L°“(IR.") (~) 
W;q^(M"). Then each Ui must be constant provided n < 4. 

Proof. Multiply both sides of (I1.15|) with RCfflui — ||wi||L~(R»)] and use assumptions to get 

(5.1) - - ||Mi||L~(R»)] div($'(|Vuip)VMi) < 0 in K". 

Applying integration by parts, for each f = 1, • • • , m we obtain 

(5.2) f $'(|Vw,n|Vw,pfiC'< 2 / $'(|Vu,^|Vu,||VRC.|[l|w^llL~(R") 

Jbr Jbr 


From the Cauchy-Schwarz inequality, for any R > 1, we get 


(5.3) 


TTl p 

Y, / < cR^-\ 

i=i Jbr 


Since u is a stable solution of (11.151) there exists a sequence of functions (j) = (c^i)^]^ that each (fi does not 
change sign. Similar to the proof of Theorem 15.11 set ifi := Vui-rj where ry = (r/', 0) £ x {0} and define 
CTi := Lemma Wff \ implies that a = (CTi)™^ satisfies (14.11) . Note that = 'f’i < iVwip. From this and 
the fact that one of conditions (A) or (B) holds there exists a constant M that is independent from R such 
that 

m p m p 

Y f>!a^A{Vu,)RT,- rT, < Y \^u,\^A{Vu,)rT,-rT, 

A — t J Br A—^ J B2R 


(5.4) 


771 p 

< mY / d>(|Vw.p), 

d B 2 R 


here we have used the fact that ||Rrj(x)||j;,cx,(Rn), ||Vui||j;,cx,(RTi) < C for some C that is independent from 
R. Note that due to the general assumption 2s4>"(s) -I- <l>(s) > 0 when s > 0 and $(0) = 0 we have 
0 < $(s) < 2$'(s)s for any s > 0. This implies that 


(5.5) 


0< /" $(|Vu*|2) < 2 / $'(|Vm.|2)|Vm. 

JBr JBr 
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From (I5.5L (15.41) and (15.31) we get 

m ^ 

(5.6) V/ nT,<CR^-^. 

i=l Jbh 

We now apply Proposition 14.II for hij = djHi{u)(l)i(pj and the identity function / to conclude that each ai is 
constant when n < 4. This implies that each m is a one-dimensional solution of p.l5p . Finally ()5.6I) implies 
that each m must be constant. □ 


In the absence of stab ility condition, there are various Liouvile theorems for solutions of (I1.15|) . at least 
for the case of m = 1, in 4^ [481 and references therein. For the rest of this section, we mainly focus on the 
p-Laplacian operator that is when $(s) = |s 2 and radial solutions of (I1.15p . For this operator, we provide 
an optimal Liouville theorem for radial stable solutions. The critical dimension is u = -\-p that is much 

higher than n = 4 given in Theorem 15.11 for not necessarily radial solutions. This implies that Theorem 15.II 
does not seem to be optimal. 

Applying the definition of the p-Laplacian operator for radial functions in dimension n, (I1.15|) reads 


(5.7) 




Suppose that u = (u^)™ i is a radial stable solution of (EZD then in the light of (11.271) and (11.261) we have 
[ |Vu,r2(V</)„VC.) + ip-2) [ |Vu,r4(Vu„V(/.,)(Vu„VC,) 

djHi{u)(l)jCi. 


= E 

i=i‘ 


R" 


We now provide the stability inequality for solutions of (11.151) with the p-Laplacian operator. This is a 
particular case of Lemma 13.11 

Lemma 5.1. Let u denote a stable solution of U.15\} . Then 

m p _ m p 

(5.8) ^ ^d,H,{u)d.H,{u)QQ <{p-l)Y. |Vu.r2|VGp, 

i,j—l i—1 ^ 

for any C = (COT w/iere Q G L°°(R”) D 1T^’^(IR") with compact support and 1 <i <m. 

For radial solutions stability inequality is of the following form. 

Lemma 5.2. Suppose that u is a radial stable solution of il.l5\) . Then 

m „ Ip,: IN m „ 

(5.9) (^-1)E/ <{p-l)Y^ / uf{\x\)\V(j){x)\^dx 




(5.10) 


E / - \j9jH,{u)diHj{u)] Ui{\x\)uj{\x\)(j)'^{x)dx, 




for all (j) G D VF^’^(R.") with compact support. 

Proof. Suppose that u = (ui)™! is a radial solutions of (11.151) that is 


(5.11) 


77, — 1 , 


/'IP-2,/ 


u'-(p-i)Krv/ = ij.(u), 


for 0 < r < 1 and i = 1, • • • ,m. Multiplying the i*^ equation of (15.111) with {u[(jPr'^ ^)' for (f G L°“(]R"') D 
pj/i, 2 (]^n) compact support and performing integration by parts we obtain 

(5.12) f ^djH,{u)u'= - f (u'(/iV"-1) 

Jr+ dR+ V ^ / 


+(!'- 1) / 

Jr+ 
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for all 0 < r < 1 and i = 1, • • • ,m. In addition, straightforward calculations show that 
(5.13) + 


(5.14) 




|u;r -u;-I = ‘^u- + {p-l)- — -\ui\P ^u''. 


Substituting (j5.13|) and (15.141) in (15.121) we get 


(5.15) 


/r+ 


djHi{u)u'ju'i(j)^r 


n —1 


n — 1 


'|P ^2 n -1 




/R+ 


u'i\P(j) r 


+ip-i) f iu'r-2u"(u'</)2)v-i 

Jr+ 


Taking sum on the index i, we get 


(5.16) 


djH,{u)u'jUi(j)^ = {n-l)J 2 






2 = 1 ' 


m p 

+(p-i)E/ 

i=l 

We now apply the stability inequality (15.811 where (j)i is replaced by u^cj) for a test function (p. Therefore, 

m „ _ m „ 

(5.17) / JdjH,{u)d^Hj{u)u',u'j(l)‘^ < (p- 1)E / ■ 

Expanding the integrand of the right-hand side we get 

(p - i)\uT-Wy{<4>)? = {p-i)\u[r^ • Vwf 

= (p - i)i<r-2 (i^x'nv(/.|2+v(<^2<) • v<). 

From this, (I5.12|l and (15.1711 we get the desired result. 

□ 


Now, we are ready to classify radial stable solutions of (E3. 


Theorem 5.2. Suppose that p,m > 1 and u is a radial stable solution of symmetric system il.l5\) where 
Hi G C'^(R'") whenever djHi{u) > 0 for all i,j = I,-- - ,m. Then, there exist positive constants vq and 
Cn,m,p such that 


(5.18) 


2=1 


\ui{r)\ > C, 


n,m,p 



ifriT^^+P, 

ifn=^+p, 


where r > rp and Cn,m,p is independent from r. In addition, assuming that each Ui is bounded for 1 < i < m, 
then n > -I- p and there is a constant Cn,m,p such that 


(5.19) 


2=1 


/ \ oo I \ h ("p+2—n+2^/ p_ J 1 

'^iy) y^i I ^ Cn,m,p'^ ^ ^ V J 


where r > 1 and := limj>_ 5 .oo Uii^) for each i. 


Proof. Let u = be a radial stable solution of symmetric system From Lemma [5^ the stability 

inequality becomes 

m „ Ip,: I, m . 

(5.20) (’^-1)E/ L P f{x)dx < (p- 1)E / ul’{\x\)\y(j){x)\^dx, 

i=i 1^1 i=i 

for all 4> € L°°(M") (h 1F^’^(K") with compact support. Note that the nonlinearity H = {Hi)”!^ does not 
appear in (15. 201) . The methods and idea that we apply in this proof are strongly motivated by the ones used 
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in 14-1^ for the case of a scalar equation, that is when m = 1, and in 2^ 3^ for the case of system 

of equations that is when m > 2. Test this inequality on the following test function (p G nL°“(R+) 


1 , 


Ht) ■= ' 


0 , 




■/. 


dz 


R 

t k'r(z)’ 


if 0 < t < 1; 


if 1 < t < r; 


if r < f < i?; 
ii R<t, 


for any 1 < r < i?. Straightforward calculations show that for the given test function the left-hand side 
of the stability inequality (15.201) has the following lower bound, 


nf IIL 

(5.21) (n-l) '^\u'f{t)e-^dt + {n-l) 

•^0 i=l -^1 i=l 


f=T+”“^ 


dt. 


Similarly we can simplify the right-hand side of the stability inequality using the fact that 

f 0, if 0 < t < 1; 

, if 1 < t < r; 


^'(t) = 


n—1 ^ Y ~ 

p-1 


'-fR — - 35 -if r < t < i?; 


.0, if i? < t. 

Substituting this in (15.201) . the right-hand side of the inequality would be equivalent to 

pr I - — 2 a 


(5.22) 


(n- 1) / t 




'^\u'f{t)dt+^ 

2=1 Jl 


/ Tt-1 

P-1 


dz 


Collecting (|5.21l) and (j5.22|) . in the light of (|5.20|) . we get 

ds 


(5.23) 


< Cr,. 


where the constant Cn,m,p is independent from r, R and it is given as 

p-1 


r ^^-1 Y,T=i l“iP(2) 


p-i VI < r < i?. 


Cn 


- 1 ) fo EZi ^dt 


Applying the Holder’s inequality we obtain 


rR 


(5.24) 


ds 


_ ^ {ET=MfZ))^ 

Ir SP+T Jr s^+T (EZi |m'|^(s))p+^ 


-ds 


< 


ji=l 

ds 


p+i / rR / 


s--^ET=M\Zs) 


/ rt / \ 


From (|5.23p we get 
(5.25) 


ds p ^ 

a-1 — ^n,m,p ' 


• r SP+i 


^ / |u'(s)Ms 

J 


Performing straightforward computation for the integral in the left-hand side of (|5.25l) and taking R = 2r, 
for any n > 2, one can get 


(5.26) 


\ui{2r) - Ui{r)\ > Cn,m,pr 


L[p+2-n+2j^] 
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Note that each m is bounded. Therefore, (15.261) implies 

m m oo 

|M,(2'=r) - Ui(2'=-V)| 

z—1 i—1 k—1 

oo , . -- 

> C^(2'=-V)Kp+2-"+2^). 

k^l 

This proves the second part of the theorem that is (15.191) and n > + p. To prove the first part of the 

without loss of generality, we assume that 2 : 

1 < ri < 2. Therefore, 

mm m 


theorem that is (|5.18p . without loss of generality, we assume that 2 < n < -\-p. Define r = 2^ ^ri where 


Z=1 


2=1 


2=1 


k-1 


i=i j=i 


2=1 


k-1 


-(p 


2 = 1 ^’ = 1 2 = 1 

This proves the pointwise bound (15.181) when 1 < n < +p. Finally, when we have dimension n = +p, 

from the above inequality, we can prove 


X l^*(^)l ^ C'n.m(fc - 1) - X 


2=1 


2 = 1 


The fact that k — 1 = log 2 ^ completes the proof. 


□ 
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